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Expecting large electric vehicle (EV) usage in the future due to 
environmental issues, state subsidies, and incentives, the impact of EV 
charging on the power grid is required to be closely analyzed and studied for 
power quality, stability, and planning of infrastructure. When a large number 
of energy storage batteries are connected to the grid as a capacitive load the 
power factor of the power grid is inevitably reduced, causing power losses 
and voltage instability. In this work large-scale 18K EV charging model is 
implemented on IEEE 33 network. Optimization methods are described to 
search for the location of nodes that are affected most due to EV charging in 
terms of power losses and voltage instability of the network. Followed by 
optimized reactive power injection magnitude and time duration of reactive 
power at the identified nodes. It is shown that power losses are reduced and 
voltage stability is improved in the grid, which also complements the 
reduction in EV charging cost. The result will be useful for EV charging 
stations infrastructure planning, grid stabilization, and reducing EV charging 
costs. 
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Indices: 


m EV index 

n network node index 
t timestamp index 

j injection index 


Variable Parameters: 


AT Nodes reactive power injection duration (hr) 


N Number of network nodes (n) 
M Number of EVs per Node (m) 


S Number of nodes for Injection(s) 


E" EV charge required (kWh) 


tst Start charging time variable (hr) 
Um Charging time variable EV (hr) 
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Pavg Average EV power constant (kW) 

P” EV active power (kW) 

PeLoss Active power loss due to EV (kW) 
P!4 Other loads on network constant (kW) 
S” App.arent EV power (kW) 

Q” EV reactive power (kVar) 

o*-ss Reactive power loss due to EV (kVar) 
Vz Network node voltage p.u. 

V” EV voltage p.u. 

Pxtm EV Power charging variable (kW) 
Antm EV charging cost constant 

n“! Optimized selected nodes 

P” Maximum charging power rate for EV (kW) 


Optimization Subscripts: 


suncoor Uncoordinated optimization 


“coor coordinated optimization 


“P Reactive power injection optimization 


‘inj, AT” Reactive power injection optimization with time interval 


‘®© Optimized values 


‘Load? Optimized values 


Optimization Parameters: 

‘’¢ Below Replaces above Optimization Subscripts 

E? Network charge (kWh) 

P’ Network active power due to injection (kW) 

Q ’ Network reactive power due to injection 

P?"5s Network active power loss due to injection (kW) 

Q *"°5s Network reactive power loss due to injection (kVar) 
V ’ Network voltage for duration AT p.u. 

P ?’6ss Network active power loss for duration AT (kW) 

OQ’ 5s Network reactive power loss for duration AT (kVar) 

f; in Optimization objective function 

{P,Q} Set of active and reactive power 


1. INTRODUCTION 

Reactive power compensation has been studied for many years with diversified research fields. The 
onset of renewable power sources and utilization brings new challenges for grid stabilization. As such EV 
load on the grid is expected to increase in the future which will aggravate electrical power supply 
stabilization and put forward new challenges to systems security, economics operation, and energy 
management. EV charging is subjective to the load on the bus network. EV charging stations and aggregators 
locations are selected on basis of lightly loaded nodes in the network buses. 

On-load tap changer (OLTC) of the main transformer, switching of the substation, feeder capacitor 
banks [1], distributed generators (DG), and adjustable loads are integrated into distribution network [2] to 
overcome losses in network buses. The uncertainty of the large number of DG may lead to over voltage or 
under voltage, increasing the volatility of the voltage, as a result, the traditional distribution network is 
gradually changing over to the active distribution network (ADN). Kotenev et al. [3] present a mathematical 
model to control reactive voltage on a node by using the synchronous motor excitation control. Meng et al. 
[4] present a key node determination method to control reactive power and voltage optimization by 
compensation devices installed at the key nodes. Zechun and Mingming [5] choose the compensation nodes 
based on the sensitivity analysis. Whereas, in [6], the reactive power margins method is proposed to 
determine the nodes where reactive power compensators are installed. Distribution feeder configuration 
(DFR) is another way that reduces power losses and manages power stability in a power distribution system, 
by sectionalizing the network by closed switches and open tie loop switches. Singh and Tiwari [7] present a 
two-stage framework is proposed for real and reactive power management and DFR. Uncoordinated and 
random charging of EVs increases peak load, losses, and voltage limit violations in the distribution system 
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[8], [9]. The uncoordinated adverse EV charging effects can be removed by coordinated scheduling of EVs in 
which an aggregator coordinates their charging/discharging with the distribution system operator (DSO). The 
aggregator act as a business entity to schedule charging/discharging of EVs and address both interests of EV 
owners and DSO in terms of economic incentives and efficient system operation of power losses and voltage 
stabilization, respectively [10]-[14]. 

EV reactive voltage stabilization with grid injection problem is usually solved in phases. The 
demand decisions of each EV are managed by regional aggregators. The objective of each regional 
aggregator is to minimize net cost and load variations by solving an optimization problem to obtain an EV 
charging schedule. Based on the real power schedule of EVs, obtained from the first phase the regional 
aggregators determine the maximum aggregated reactive power that can be injected/absorbed by all EVs 
present at each node. This information associated with maximum reactive power injection/absorption of each 
node is sent to the DSO by each regional aggregator as shown in Figure 1. Thereafter, in the second stage, 
DSO performs simultaneous reactive power scheduling with power flow calculations to minimize losses in 
the distribution system. Wang et al. [15] presents coordinated EV charging with reactive power support to 
distribution grids at each network node. Shafiee et al. [16] investigate the impact of plug-in hybrid electric 
vehicles (HEV) on power distribution systems but didn’t offer a solution for grid power losses. Modern 
inverters have the capability of providing reactive power support to improve voltage profiles and to minimize 
power, several authors have addressed this topic using the bi-directional EV battery charger enabled with the 
vehicle-to-grid (V2G) concept such as by Yong et al. [17]. The EV chargers also have the capability of 
allowing reactive power flow to the grid by utilizing the direct current (DC)-link capacitor. 


‘Distribution Gid Tg 


Figure 1. IEEE 33 node power network [15] connected with EV aggregator chargers and distributed grid 
control centre/DSO 


The novel approach adopted in this work is to first search the network nodes which are most 
affected in the grid due to incoherent EV charging and then optimize the reactive power injection magnitude 
and time duration. Whereas all the research work mentioned above; studies the EV charging modeling and its 
impact on the grid and offers a solution to reduce the grid power losses by reactive power injection on all 
nodes irrespective of injection time duration. A large-scale 18K EV charging model is simulated in this work 
with driving pattern, EV charging intervals distributions, and EV charging cost. IEEE 33 node radial network 
[18], [19] is used for the study. This paper aims to regulate the grid voltage, and a new approach is 
considered here to minimize the EV charging voltage rather than the traditional minimization of the grid 
power losses approach. From an optimization point of view, the minimization decision criteria are based on 
the mean square error (MSE) of grid voltage while the input control variable is reactive injection power. 

A sequence of numerical computations adopted in this work is described in sections 2-4. EV 
charging parameters are defined in section 2 and Table 1. EV user charging requirement is modeled with 
Weibull distribution in section 2.1. The daily demand profile and cost of electricity used in this study are 
discussed in section 2.2. In section 3.1 uncoordinated EV charging load without daily load, the profile is 
evaluated with B/F sweep power flow analysis. In section 3.2 the coordinated EV charging load and charging 
cost is optimized with load profile and charging constraints using mixed integer linear programming (MILP). 
In section 4 most affected power losses and voltage instability nodes are searched, followed by the 
optimization of reactive power injection magnitude required for its restoration. Finally, in section 5 the 
reactive power injection time duration is determined at the located nodes. 
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2. EV PARAMETERS AND CHARGING 

IEEE 33 network [19] data in appendix 1 is used for simulation studies in this work. A large number 
of EVs=18000 are used in the simulation for maximum EV impact on the grid. Each node can act as an EV 
aggregator and can carry max M = 563/node EVs with different state of charge (SoC) and random 
charging interval. Simulation studies are carried out for 24 hrs over 15 min time interval (96 time stamps). 


2.1. EV charging distribution 

Due to different brands of EV cars with battery energy capacities Ea ranges between 4kWh - 
50kWh and considering local driving pattern the energy requirement probability is represented by the 
Weibull density [18] curve defined in (1). 


ij 


a 


(1) 


re b-1 

req) — 2 (e Z 
fE )=2.( a ) e 
Where choosing a= 15, b= 1.4, c= 2, the probability curve is given in Figure 2. Nominal average EV 
charging is considered as Pavg = 7.2kW. EV start charging times t,, are randomly generated based on local 
EV user charging behavior pattern shown in Figure 3 is based on 20% user charging between 07h00 - 10h30, 
40% between 16h00-21h00, and the rest is evenly distributed over the day. The time interval U,, for mth EV 


is then random, both in start time t,, and energy required ie and is defined by (2). 


Um = (tse; stop time] 


Enea 
Um = is tst + | (2) 
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Figure 2. EV charging energy requirement is modeled with Weibull density probability 
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Figure 3. EV charging start time tst histogram 


2.2. EV charging and network parameters 
To simulate the EV charging process parameters in Table 1 are selected. General equatorial 
condition daily demand pattern and cost of electricity data used in this study are shown in Figure 4. The 
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backward/forward sweep method [20] for power flow analysis is used for IEEE 33 radial distribution system 


shown in Figure 1. Baseload active and reactive power (P, Q)4°% , (P, Gy een and voltage (Vz) p.u. of 
TEEE 33 radial system [19], is given in appendix 1. 


Table 1. Network and EV parameters 


Network PQ nodes, N 32 

EV per node, M 563/node 
Max EV charging Rate Pon 9.4 kW 
Average EV Charge Pavg 7.2 kW 
Slack Gen PV node 1 5 MW 
EV charger PF 0.74-0.98 
Timestamp At=/5min 96 min-step 
Total EVs (M . N) 18,000 
Base MVA 1000 
Power factor, Phase Angle 0.74, -42.30 


Daily Load Profile 
Electricity cost 


Figure 4. Daily load profile and cost of electricity 


3. EV CHARGING 

The high volume of EVs charging with the grid affects the grid voltage and power factor. The 
capacitive load of EV grid charging impacts the grid [21]. EV aggregator is expected to meet the demand of 
battery charging at an economical cost. EV power charging is modeled and simulated as a function of three 
variables EV number (m), Network node (n), and time (f). 


3.1. Uncoordinated EV charging 

For uncoordinated EV charging Entm (kWh) is obtained from the Weibull distribution in (1) and 
charging interval Uņ,tm is evaluated from (2). The EV charging power P;% required for all M=563/node on 
32 nodes (N) for 15min time step t (1:96) is given in (3). 


Zejo tne 
Pot = Lm? Go (3) 
n,t,m 


EV charger power factor [22] is varied between pf = [0.98: 0.74 ], pf=0.74 corresponds to a —42.3° angle 
between the voltage and current phasors. The maximum reactive voltage over time is shown in (4). 


uncoor — ev 
Pi = max[Pre] (4a) 
uncoor — [Patter 
sr = SE w 
Quncoor = [sper Ma Puncoor? (4c) 
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3.2. Backward/forward sweep power flow 

The backward/forward sweep method [16], [19] for power flow analysis is used for power flow 
analysis. Active power P#"°°°Tand reactive power Q7”°°°” flows in a branch between node ‘n’ to node 
‘n+l’, with branch resistance R, and impedance X,, given in appendix 1. B/F sweep method calculates 
backward from (n+/) node the last node (n) is given in (5) and (6). 


uncoor! punçoori? + qunçoorı? 
uncoor — n n 
Pr = Pasi tin — AA (5) 


Z 
Vh+1,t 


2 2 
uncoor! uncoorl 
Pn+t +Qn+i 


uncoor — juncoor: 
Qn Zo Oni + Xn 


(6) 


Viteat 
Where the update values are, Pup {O0"” = Pungor + [phodd + PSY] 


aoe = Quran 4.) Load y + Qe v] 


Prod and Q}95f are IEEE 33 network loads in appendix 1, P$% and Q$% are the EV loads that are 
connected at node ‘n+1’. P¥peoor | Quncoor and V,2., are effective real, reactive power flows in network 
branches due to loads. The network node voltage is given in (7), 


Vater SVN TE Gg) (7) 
Power Losses are shown in (8) and (9) 


puncoor.Loss = I2 "a Rn (8) 


n 
greer kass = I *jXn (9) 


the voltage magnitude V""°°°" and the angle at each node are calculated in a forward direction and are used 
recursively in a forward direction to find the voltage and angle respectively of all nodes of radial distribution. 


3.3. Coordinated EV charging with constraints and cost-MILP 


The interest for the EV aggregator is to minimize the charging cost. Considering Px, as decision 
variable for EV charging power, and A, + as charging cost. To optimize the EV charging load according to 
the charging cost and the load profile, the objective function f°9°" minimize the product of EV charging 


power and cost of charging as (10), 


coor _ 

min =.) par Px tım Antm (10) 
O<ts9 
demem 


subjected to constraint, 


Pxem < P™ (Max charging Limits) (11) 
Prem < [P$}m + PL4] (EV + Base Load) (12) 
Pin Unie = Ben (Energ Balance) (13) 
Patm 20 (14) 


constraints in (11), (14) are maximum charging power limit, optimized charging should be less than or equal 
to total inclusive of node load, optimized charging energy should be equal to required EV charging energy 
and positive charging power respectively. 

Due to a large-scale optimization problem, IBM ILOG CPLEX optimization linear programming 
function ‘cplexip’ is used in MATLAB 2019 environment. For grid analysis, the maximum load value is 
considered at the node. Total nodes (n) load at t + At, where At =/5min, is given in (15), 
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* 


coor — Y'‘M=500 pev* . ev* _ px 
Prt = es Prim ne Pritam = Fntm (15) 


where Paes is the coordinated optimized active power and QEY m can be found out similarly to (4). The 
maximum EV charging power at each node is defined in (16), 


E= maa PEA a] (16) 
N€32 


the Backward/Forward sweep power flow analysis is used to evaluate the network nodes parameters, 


coor — ev* 
Pree" = max[Prie | (17a) 
geoor — PA” (17b) 
i Pin 


Qer = | scoor? — pcoor? (17c) 


2 2 
coor! coor! 
P, + Qn41 


ppor = Profe Ro (18) 
n+1 
coor — pcoorr peger'? + ogoori? 
ot Oe (19) 
where, 
peoort = pooor + ptega 
A A i 
and, 
VE = Yeoor T(r, + jXn) (20) 
preoa Loss = I Ra (21) 
QRS = TR IX (22) 


4. REACTIVE POWER INJECTION OPTIMIZATION-GA 

Power losses and voltage stability can be improved with reactive power injection in (6). The 
optimized EV power a (S) nodes are searched here by GA to reduce power losses and stabilize the 
voltage. The reactive power may come from shunt capacitor bank, EV battery DC-link, or offload 
synchronous motor condenser. It is necessary to introduce injection power magnitude variable term Q[j,’] as 
a function of node location index j,,” in (24) for mixed-integer GA optimizer as, 


Pp’ = max[Pey" 
i) = max PEE 


inj 


oii as Pn 
n Pin 
ae in 2 in 2 
inj _ inj inj 
n 7 Sh = Pa 
inj? ninj!” 
inj inj’ pnt +Q 
mj _ plnj n+1 n+1 
Pn Sed + Ry v2 (23) 
n+1 
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injr? | „inj? 


Y= Olin] + Oh + Xp EE pea ea 
where, ~. 7 
Pia =R” Pee 
i 2 inj + Qload 
and, 
pile) (25a) 
WA = 1 pu. 
Vang = Va — In” (Rn + Xn) (25b) 
pines 2 ri oR. (26) 
inj Loss _ co JX (27) 


the reactive power injection variable Q[j,] is defined as the customized function for ‘Population Generation’. 
The index s = 1-:S, are number nodes which optimizer will be searched to minimize the objective function 
(29). Selecting S=3 sets the remaining (32-S) nodes are set to zero. Q[j,] denotes the reactive power injection 
magnitude variable, also searched by the mixed integer genetic algorithm (GA) optimizer. 


js = randperm{ 0, 0, 0, Imad, 0, O, Ln p- --0,..-, Inci7s-+s0]ix33 (28a) 


0(1), 0(2) ,0(3), Q[a'] ,0(5),0(6), QUi7'] pose 
Ql] = Q ,0(12), cig OL pat iy 033 1x33 (28b) 


The reactive power injection Q[j,] is assigned a range of [0 300]KVar to minimize the objective function in 
(28)b. MATLAB command ‘randperm’ is a random and combinatorial set of cyclic permutations for GA 
optimization. Minimum voltage objective function is formed with two voltages, the constraint Vz IEEE base 


voltage p.u and variable term ed (25(b)) to represent the EVs charging voltage which is a function of I a 
P.” and Q,,” in (23-27). 


inj 


pin} | [ENV] + Max(vir’)’ (29) 


4.1. Reactive voltage injection time duration interval-GA 

The optimized EV power pev’ for S=3 nodes n°” = [14; 17; 18]3,1 selected in section 4 injects 
the reactive power over a 24 hour period which may not necessary for dynamic EV loads. To find optimal 
injection time interval (AT) hour and reduce injection charges at nodes nse = [14; 17; 18]3,1, the 
optimization in section 4.0 is repeated with timestamp (t) variable and B/F sweep power flow analysis with 
time function is, 


inj,AT _ pev* 

Pit —~ Ant 
inj,AT 

sini AT _ Pht 
nt Pint 


inj,AT 


on Ame 
inj,AT 
nt S 


ini Ame 
inj,AT 
nt =P 


n,t 
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inj ATI? inj ATI? 
inj,AT _ pinj,AT' Patist  +2ntit 
Pat i Patist + Rn v2 (30) 
n+1,t 
inj ATI? | inj, ATI? 
inj AT _.. inj, AT' ntit  +On+it 
nt N gyset + Onsite + Xp Vast (31) 


where the updates with time function (t) are, 


inj,AT' _ pinj AT Load 
P = Pir + Prit 


nt+1,t 
inj,ATr _ ,inj,AT Load 
n4+1t — ‘nt + Qnit 


using GA MILP optimization with customized random ‘PopulationGeneration’ hourly timestamp is 
randomly searched over 24 hours in three nodes nse; defined as, 


0, 0, 13 =e 024 242 
injtn = rand |1, 02 03 + O24] }224 
0, 0, 03 base 124 102 
and, 
inj,AT inj,AT . 
Vii = Var — int: (Rnt + jXne) G2) 
pee = poer „Rat (33) 
genes = par JXnt (34) 


the Objective function is defined as before in (29) with an additional time interval (t) variable to determine 


the reactive power injection Q[ns°"]*. The minimum time interval for the selected nodes nS! = 
[14; 17; 18] yr. 
inj inf AT\2 
fain = Jx V7 + Malye) (35) 
5. RESULTS 


Coordinated power EV charging P{°°" and 1hr injection EV charging power P,” AT are shown in 
Figure 5, with higher power utilization at lower daily demand and electricity rates. Uncoordinated [23] EV 
charging the penetration to the grid starts immediately when the user arrives at home and plug-in. Most EV 
users arrive at home at the same period of peak demand creating a large load, power losses increased 10 
times as shown in Figures 6-7, and voltage instability increases in Figure 8. Coordinated charging with the 


objective function fpin and constraints of variable P%tm in (11) and (14) reduces the network power losses 


pooenboss QcoorLoss in Figures 6-7 and stabilize the voltage V,°°°" in Figure 7. GA optimization with the 
objective function je for power losses P” eng ee SS and voltage drop y” which randomly search js 


nodes for 24hr reactive power injection results in selected nodes jž- [n$®%]sxı are shown in Table 2. The GA 


aes a 7 5 2 inj,AT inj,AT,Loss inj,AT,Loss 
optimization objective function f% pe 7 


‘min for Ihr reactive power injection results in 5 Qn ; 
inj,AT* 
Vn 


further reduces the active and reactive power in Figures 6-7. 


Table 2. Three node injection 


Inj. Voltage Inj. Hr 

Q[ ns" }* Jenset 
242 1 
224 1 
102 18 
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Figure 5. Daily load profile and cost of electricity. Shows the excess EV charging power required in 1-7hr 


Active Power Losses KW 


Figure 6. Active power losses 


Reactive Power Losses Kvar 


Figure 7. Reactive power losses QŁ944, Qf 


It is observed 
reactive power losses 


power injection {P,Q} 


with coordinated and 1hr reactive power injection 


—+*— 24hr. Inj. Power Losses 
——*— Uncoor. Pow. Losses 
—— 1thr.Inj Pow. Losses 
—*— Coor. Cost Pow. Losses 
IEEE Base Pow. Losses 


puncoor,Loss při oss 


Load 
Pr sin sin 


coor,Loss 
+1” Pa 


inj,AT,Loss 
and P” 


—+*— 24hr Inj. Paver Losses 
—@—-  Uncoor. EV Pow. Losses 
—s— thr.Inj Pow. Losses 
—=—— Coor. Cost EV Pow. Losses 
IEEE Base Pow. Losses 


20 


35 
Node 


coor,Loss uncoor,Loss 


inj,Loss 
Qn > Qn an 


inj,AT,Loss 
d Qn 


that 1hr injection traces closely matched with IEEE Base network traces. Active and 
of IEEE baseload {P, Q};077"°% {P,Q} EP, QRL, 24hr reaction 


inj,Loss ý PE : inj,AT,Loss . š 
ma and 1hr reaction power injection {P, Q}? are shown in Figures 6-7. 
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IEEE network p.u. voltage Vz, V09" , Vreo or Wie and Vee AT are shown in Figure 8. 1hr voltage trace is 


the most stable voltage matching with IEEE network voltage Vz. Table 3 summarizes the comparative results 
of the techniques used in this research work. 
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; coor yuncoor inj inj,AT 
Figure 8. Nodes voltages p.u. VL, Vy VY Vie and Ve 


Table 3. Results of techniques used 


Analysis Sum of Active Power Losses Sum of Reactive Power Losses Voltage MSE 
į (kW) (kVar) 
Uncoordinated EV Charging 477.72 329.06 7.6568 
Coordinated EV Charging 69.51 46.62 0.7993 
24 hr. Injection EV Charging 56.04 37.53 1.4417 
1 hr. Injection EV Charging 46.02 30.77 0.3952 
IEEE Base load without EV 28.15 18.79 0 


Despite stabilizing the grid voltage and improving power losses at the most affected nodes. The 
optimized power injection magnitude has the drawback of requiring all the nodes to be equipped to inject 
reactive voltage back V2G into the grid. This technique is efficient for known injection nodes at the hot spot 
(commercial outlets, offices, restaurants) injection nodes [24], [25] and equipped with reactive voltage 
injection devices. 


6. CONCLUSION 

A large number of EVs charging with the grid affects the grid voltage and power factor. EV grid 
charging impacts the grid with capacitive load. In this paper, the grid voltage recovery is demonstrated for 
18K EVs charging load due to load profile capacity. It is shown that optimized reactive power injection 
magnitude and time interval at the most affected nodes stabilize the grid voltage and improve power losses. 
Changing EV's user driving pattern distribution, cost (A), and local load profile, will help to locate the 
installation of EV charging stations and the aggregators to inject reactive power for power stability. EV 
charging analysis presented in this work can be extended and compared with distribution Feeder 
reconfiguration (DFR) of the IEEE network to study the grid stabilization problem. Further, the number of 
EVs can be increased and analyzed with higher profile loads. 


APPENDIX 
IEEE 33 BUS DATA 
; ; a ; Load at receiving end Bus Lewd.L Load L 
Ras Se a i a T Real Power Reactive Power He ae es as 
umber us us (kW) (kVAr) (kW) (kW) 
1 1 2 0.0922 0.0477 100.0 60.0 16.80 8.57 
2 2 3 0.4930 0.2511 90.0 40.0 71.39 36.36 
3 3 4 0.3660 0.1864 120.0 80.0 27.67 14.09 
4 4 5 0.3811 0.1941 60.0 30.0 26.04 13.26 
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IEEE 33 BUS DATA (cont.) 
; : on ; Load at Receiving End Bus madtoss AL oad Lbs 
ie ae ae aes S ae sa aad Real Power Reactive Power en a 
¢ : (kW) (kVAr) 
5 5 6 0.8190 0.7070 60.0 20.0 53.30 46.01 
6 6 7 0.1872 0.6188 200.0 100.0 2.68 8.85 
7 J 8 1.7114 1.2351 200.0 100.0 6.79 2.24 
8 8 9 1.0300 0.7400 60.0 20.0 5.89 4.23 
9 9 10 1.0400 0.7400 60.0 20.0 5.02 3.56 
10 10 11 0.1966 0.0650 45.0 30.0 0.78 0.26 
11 11 12 0.3744 0.1238 60.0 35.0 1.24 0.41 
12 12 13 1.4680 1.1550 60.0 35.0 3.77 2.97 
13 13 14 0.5416 0.7129 120.0 80.0 1.03 1.36 
14 14 15 0.5910 0.5260 60.0 10.0 0.51 0.45 
15 15 16 0.7463 0.5450 60.0 20.0 0.40 0.29 
16 16 17 1.2890 1.7210 60.0 20.0 0.36 0.48 
17 17 18 0.7320 0.5740 90.0 40.0 0.08 0.06 
18 18 19 0.1640 0.1565 90.0 40.0 0.21 0.20 
19 19 20 1.5042 1.3554 90.0 40.0 1.11 1.00 
20 20 21 0.4095 0.4784 90.0 40.0 0.13 0.16 
21 21 22 0.7089 0.9373 90.0 40.0 0.06 0.08 
22 22 23 0.4512 0.3083 90.0 50.0 4.30 2.94 
23 23 24 0.8980 0.7091 420.0 200.0 6.96 5.49 
24 24 25 0.8960 0.7011 420.0 200.0 1.74 1.36 
25 25 26 0.2030 0.1034 60.0 25.0 3.65 1.86 
26 26 27 0.2842 0.1447 60.0 25.0 4.68 2.38 
27 27 28 1.0590 0.9337 60.0 20.0 15.90 14.02 
28 28 29 0.8042 0.7006 120.0 70.0 11.03 9.61 
29 29 30 0.5075 0.2585 200.0 600.0 5.49 2.80 
30 30 31 0.9744 0.9630 150.0 70.0 2:25, 2.23 
31 31 32 0.3105 0.3619 210.0 100.0 0.30 0.35 
32 32 33 0.3410 0.5302 60.0 40.0 0.02 0.03 
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